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- - Abstract 
OO 

\ Light in a dielectric medium moves slower than in vacuum. The corre- 

CN ' spending electromagnetic field equations are then no longer invariant under 

Q^. ordinary Lorentz transformations, but only under such transformations cor- 

\ responding to this reduced velocity. Based on this physical symmetry, an 

effective theory for low-energy electromagnetic phenomena in dielectrics is 
constructed. It has none of the problems of the old formulations of Minkowski 
and Abraham. Dispersion in the optical regime and the Kerr effect arise in 
^ ^ a natural way from higher order interaction terms. The effective field theory 

Q_i' is quantized by standard methods and quantum corrections can be calculated 

in a systematic way. Thus it relates many different classical and quantum 
^ \ optical phenomena into a unified description. 

^ . 
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Light in a isotropic, dielectric medium with refractive index n > 1 moves with a 



> 

^ ' wave velocity 1/n when setting the velocity of light in vacuum to c = 1. The 

Maxwell equations describing the electromagnetic field in matter are well known 
^ . and one would think that a consistent, theoretical description of electrodynamics 

^ ! in dielectrics would exist. But for close to a hundred years two different theories, 

O I one due Minkowski and the other to Abraham, have been used. They were both 

constructed to be consistent with the special theory of relativity in vacuum[T]. 



The main difference between the two formulations is found in the energy-momentum 
^ ■ tensors. In the Minkowski theory this tensor is not symmetric and thus has problems 

■ with the conservation of angular momentum. Abraham instead imposed such a sym- 

metry with the consequences that the momentum density of the electromagnetic field 
is reduced by the factor and a corresponding new volume force appears. When 
quantized, the photon in a medium has a non-vanishing four-momentum squared. 
Many attempts have been made to understand the underlying problems without 
any real success[2]. In standard textbooks the Abraham formulations seem to be 
preferred [3] [H, while most experiments favor the Minkowski formulation [5j. How- 
ever, for a few experiments where the system undergoes acceleration, the Abraham 
theory might be required[6]. 
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In the following it is pointed out that covariance under vacuum Lorentz transfor- 
mations is not consistent with the symmetries of the underlying Maxwell equations. 
Instead we propose an effective theory based on Lorentz invariance under trans- 
formations corresponding to the physical speed of light 1/n in the medium. We 
then avoid the old problems of Minkowski and Abraham. Since the medium itself 
is not invariant under these transformations, this new formulation is to be used in 
the rest frame of the medium. It is thus on an equal footing with other theories 
in condensed matter physics for excitations or quasiparticles with linear dispersion 
relations. It is also shown how non-linear dispersion and the Kerr effects follow from 
higher-dimension operators in the effective theory. 

The standard theory for electromagnetic fields in a medium without any free currents 
or charges is given by the ordinary Maxwell equations V x E + 9B /dt = and 
V X H - dB/dt = together with V ■ D = V ■ B = where D = and B = /xH 
for isotropic matter. This is already an effective description on large scales where 
all the underlying atomic physics describing the material content of the medium is 
parametrized by the electric permittivity e and the magnetic permeability /i. We 
will use units so that in vacuum both Eq = 1 and fiQ = 1 so that the speed of light 
c = 1/ y/eofiQ = 1 in vacuum. The electric and magnetic fields have then the same 
dimensions. 

In the following both e and ^ are taken to be constants. We thus first consider phe- 
nomena in the dielectric within a narrow range of frequencies where this assumption 
holds. The above Maxwell equations then define our field theory to lowest order. It 
is defined in the rest frame of the medium. Taking curl of the first Maxwell equation 
and using the second equation, one then obtains in the usual way the wave equation 

(^'^-V^)E(r,t) = (1) 

for the electric field and the same for the magnetic field. The electromagnetic wave 
velocity is therefore 1/n — 1/ y/ejl in the medium. 

Neither this equation nor the original Maxwell equations are invariant under the 
ordinary Lorentz group. Instead they are invariant under what we will call the ma- 
terial Lorentz group, with the vacuum light velocity replaced with 1/n. The original 
treatments of electromagnetism in matter were generally formulated covariantly so 
to maintain invariance under the vacuum Lorentz group. We do not see how that 
can be justified for phenomena on a macroscopic scale where this does not represent 
an inherent symmetry of the system. It would not be the first time in physics where 
the symmetry is different on different scales. 

It is straight forward to construct a special theory of relativity within the medium. 
One just has to replace the ordinary vacuum light velocity with the dielectric velocity 
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1/n. Taking the metric to be the same as in vacuum, a coordinate four- vector would 
then be = {t/n,x.) so that the corresponding dielectric four-momentum becomes 
= {nE,p). The equivalent of a massless photon in the ordinary sense within 
the medium should therefore have energy and momentum related by the invariant 
P'^Pfi = OT E = p/n. This is satisfied by the photon in the Minkowski theory. In 
solid state physics this is usually said to be an excitation or quasiparticle with a 
linear dispersion relation. 

A covariant formulation based on this modified Lorentz group also follows naturally. 
Combining the electric potential $ and the magnetic vector potential A into the 
four- vector potential A'^ = (n$,A), the components of the electric field vector 
E = — — dA/dt and the ones of the magnetic field B = V x A form the 
antisymmetric tensor F^^, = d^A^ — d^A^. The invariant Lagrangian density for the 
electromagnetic field within the dielectric is then 



From here one finds the standard energy-momentum tensor in the medium to be 
given by /iT^^ = F^°'Fai, + {1/ A)5^^,F'^p which is obviously symmetric. Its conserva- 
tion d^T^j^ = follows from the Maxwell equations of motions. The time component 
of this conservation law gives d£ / dt + 'V ■ N = where £ = (l/2/i)(n^E^ + B^) 
is the energy density of the field and N = E x H is the Poynting vector. This 
equation thus gives energy conservation. Similarly, the spatial components give the 
vector equation dG /(9t + V ■ T = where the three-tensor T contains the Maxwell 
stress components Tij. We therefore conclude that the vector G = D x B = n^N 
represents the momentum density of the field. 

Quantization of this free theory is straight forward. Integrating the above energy 
density, the resulting Hamiltonian gives a photon with wave number k an energy 
where = \^\/n. Similarly, the momentum density G gives it a momentum 
hk. The same energy and momentum also follow from the Minkowski theory. But 
while we would say that this particle in the medium is massless, in the Minkowski 
theory it has a negative mass squared, i.e. it is a tachyon. 

There are two obvious experimental consequences one can extract right away from 
this free theory. One is the Casimir force |7], which now can be measured with good 
precision for the geometry of two parallel plates with separation L in vacuum [8]. 
Replacing the vacuum with a dielectric the force is expected to change. This experi- 
mental situation has recently been considered by Brevik and Milton [Qj who after 
a rather long calculation find within the Minkowski theory that the Casimir force 
should be reduced by a factor n compared with the vacuum force F = — /i7r^/240L'^. 
This result now follows straight away from the above since the force is essentially 
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just the zero-point energy E = J^k^^k between the plates. The sum Eq = J2kM^\ 
represents the corresponding zero-point energy in vacuum and is independent of the 
refractive index. It is not clear what the alternative Abraham theory will give for 
this Casimir force. 

We can also make a similar, direct derivation of the energy density for black-body 
radiation in a cavity filled with a dielectric at temperature T. It follows from the 
usual expression 




with /3 = l/ksT where fc^ is the Boltzmann constant. Again with uj = k/n this 
gives in the large-volume limit a result larger than the standard Stefan-Boltzmann 
formula by a factor n^ITU]. For a typical dielectric this increase will thus be in the 
range 2-7 and its effect should be detectable provided that this lowest-order theory 
contains the dominant physics. 

Higher order corrections to these results can now be calculated within an effective 
field theory. Such theories have a long tradition in condensed matter physics and 
are now also used extensively in low- and high-energy particle physics. They are 
valid below an upper energy cut-off characteristic of the problem under investiga- 
tion and involve only the relevant field variables of lower energy. All degrees of 
freedom occurring on energy scales above this cut-off are integrated out to give 
higher dimensional couplings between the low-energy fields. Even if the theories are 
not renormalizable, quantum corrections can be calculated in a systematic way with 
divergences absorbed into the coupling constants. 

In order to be gauge invariant, higher order couplings can only involve the fields E 
and B and derivatives of them. For the sake of counting, we can assign the dimension 
+2 to each of these fields and an increase in dimension by -|-1 for every derivative 
using quantum units where h = 1. To be invariant under time- reversal and ordinary 
rotations, such new couplings must involve at least two space-time derivatives. For 
example, one possibility could be the term E ■ 9^E with 9^ = d^d^ = n^d^ — V^. 
It has dimension 6 and is also invariant under the material Lorentz group. But the 
lowest order equation of motion ([1]) is just d^E = and this term can therefore not 
contribute. Possible new terms of dimension 8 would be (E-E)^, (B-B)^, E^B^ and 
(E ■ B)^. Requiring Lorentz invariance, the first three of these non-linear couplings 
must appear in the combination (n^E^ — B^)^. All such terms represent anharmonic 
interactions involving four fields. 

Dispersion can now be obtained from dimension-6 interactions when we break the 
material Lorentz invariance down to rotational invariance. One example is VjE- ViE. 
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It is equivalent to E ■ V^E by a partial integration in the action integral where it 
appears. The similar term dtE ■ dtE involving two time derivatives is for the same 
reason equivalent to E ■ V^E when we use the equation of motion. An interaction 
like E ■ V^B is ruled out by parity invariance. 

In dielectric materials magnetic effects are negligible, and it is therefore reasonable 
to assume that all the terms involving the magnetic field, can be dropped. Including 
interaction terms up to dimension-8 operators, we then have the effective Lagrangian 

C = l{n^E^ - B^) + . V^E + A(v^e)^ + ^(E ■ E)^ (3) 

The parameter M is an energy scale which characterizes the microscopic physics we 
have integrated out. In a dielectric this is set by atomic physics and should have a 
value around 5-10 eV. It is possible to calculate or at least estimate the dimen- 
sionless coupling constants di, d2 and a from the atomic physics of the material. 
With a realistic value for the energy scale M, the Lagrangian contains only three 
independent parameters. For each material they can therefore be determined by 
three different measurements. It would then be able to predict the outcome of many 
other experiments. In practice, we expect the term oc 62 to be dominated by the 
di -term. 

Let us first consider the effect of the dispersive term oc di in the effective Lagrangian. 
It gives a contribution to the classical equation of motion for the field that is linear. 
For a plane wave of the form E oc e^e^-''-'^*) we obtain a dispersion given by u = 
k"^ /ri^ + 2diuj'^k? / eM"^ . Since the last term is assumed to be small, we can there 
set uj = k/n. With e = ra^ in a dielectric, we find the refractive index, defined 
by the phase velocity u/k = l/n^u), to be n^u) = n{l — diu^/M"^). This is of 
the same form as the phenomenological law due to Cauchy and usually written as 
n(A) = A + i?/A^[TT]. The measured values of the dispersion parameter B for most 
gases and transparent dielectrics are in the rather narrow range (2 — 15) x 10~^^ m^. 
The coupling constant di must therefore be negative to give normal dispersion. 
Comparing with our result, we find that these measured values correspond to taking 
nM = 10 eV and -di in the range 0.1 — 1.0. This is just as expected from the scale 
of the underlying atomic physics . Including the d2 -term, we would similarly have 
have obtained a correction to the above dispersion relation. 

The four-field coupling in ([3]) can describe many non-linear phenomena in dielectrics. 
It has a coupling constant a which corresponds to the third-order susceptibility x^^^ 
in the language of non-linear opticsp^. In the DC Kerr effect an external, electric 
field E affects the propagation of light. The refractive index becomes field depen- 
dent. Here it results by replacing one factor E ■ E in the interaction term with the 
square E"^ of the external field. The remaining, oscillating field then obeys a linear 
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equation of motion as above giving the refractive index n{E) = n+2aE'^ /nM"^. Mea- 
surements are often parametrized as n{E) = n + XKE"^ where A is the wavelength. 
Here K is the Kerr constant which we thus find to be given by XK = 2a/nM^. For 
isotropic glasses it has values of the order 10~^^m/V^ corresponding to a coupling 
constant a ~ 10~^ when we take nM = 10 eV. For dielectric liquids this effect is 
usually larger. 

In the AC Kerr effect the field of the wave itself affects its propagation. Using a 
semi-classical approximation, we can again replace one factor E ■ E in the interaction 
with the intensity I = eE^ of the beam. Then we find for the refraction index n(J) = 
n{l + al/ (nMY). Writing it on the form n{I) = n+n2l, we thus have n2 = a/n^M'^. 
Measurements in optical glasses [T2] give values for n2 of the order 10~^° m^/W. With 
nM = 10 eV this is consistent with a ~ 10~^ as obtained above. In fact, eliminating 
a between these two Kerr constants, we find XK = 2n^n2. This relation seems to be 
roughly satisfied by data both for optical glasses and for dielectric liquids. A similar 
relation can also be derived by simple, atomic arguments [12]. 

With this effective field theory we can now also calculate quantum corrections to 
some of the processes previously discussed. For the energy of black-body radia- 
tion in a dielectric the calculation would be very similar to what has already been 
done [13] based upon the corresponding effective Euler-Heisenberg theory for low- 
energy QED[T3]. In that theory the electron mass rrie plays the role of M here and 
the dimensionless coupling constants are powers of the fine- structure constants and 
thus very small. Since me/M k, 10^, the higher order corrections we consider in a 
material would therefore be enormously bigger than these QED effects and therefore 
in principle experimentally more accessible. But here the dominant contribution will 
come from the di dispersive term. The corresponding term does not exist in the 
Euler-Heisenberg effective Lagrangian because that theory is required to be Lorentz 
invariant. It will give a correction with size of the order (fc^T/M)^ times the lowest- 
order result. Taking M = 10 eV, this would first start to be significant at such high 
temperatures that the material has already melted. But here the lowest-order result 
itself is important to verify experimentally. 

Corrections to the Casimir force can also be obtained analogously to calculations 
previously done in vacuum [T5]. but now based on the more important di dispersive 
term. However, for the Casimir force it is possible that a surface term will also be 
present in the effective theory [T6] [17j. The dominant correction will then be of the 
order 1/LM times the vacuum force, provided that the scale of this interaction is 
of the same order of magnitude as in the bulk. It will then begin to be important 
for plate separations approaching the cut-off distance L = 10 nm. With modern 
nanotechnology this should soon be within reach of experiments . For even shorter 
separations the effective theory is not expected to be valid anymore. 
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Looking back at previous theories, one sees at least two reasons for their difficuhies. 
The first is that they were covariantly formulated based on the vacuum Lorentz 
group, a symmetry which is not present in the macroscopic Maxwell equations. In 
the present proposal the rest frame of the medium is a preferred frame. Secondly, 
the problem with quantizing them was frustrated by the effect of dispersion. Here 
dispersion is not part of the lowest-order theory, but arise as a higher order effect 
which is included using ordinary perturbation theory. 

We can also write down corresponding effective Lagrangians for electromagnetic 
fields in anisotropic materials. The rotational invariance so far assumed will then 
be replaced with the invariance under a discrete, crystal group appropriate for that 
particular material. As a consequence, the coupling constants appearing in ([3]) will 
be tensors and new couplings can arise. Needless to say, there will then be a few 
more unknown parameters, but the effective theory might still be useful. 

This work has been helped by discussions with I. Brevik, C. Burgess, T. Curtright 
and L. Mezincescu. Colleagues at the Department of Physics at UM are acknowl- 
edged for their hospitality. 
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